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In the framework of the linear sigma model, we study the time evolution of a system of classical 
' a and pion fields coupled to quarks. For this purpose we solve numerically the classical transport 

Q\ . equation for relativistic quarks coupled to the nonlinear Klein-Gordon equations for the meson fields. 

■ We examine evolution starting from variety of initial conditions corresponding to spherical droplets 

of hot quark matter, which might mimic the behaviour of a quark plasma produced in high-energy 
, nucleus-nucleus collisions. For large droplets we find a strong amplification of the pion field that 

oscillates in time. This leads to a coherent production of pions with a particular isospin and so 
would have similar observable effects to a disoriented chiral condensate which various authors have 
suggested might be a signal of the chiral phase transition. The mechanism for amplification of the 
pion field found here does not rely on this phase transition and is better thought of as a "pion laser" 
which is driven by large oscillations of the a field. 
PACS: 25.75.+r, 11.30.Rd, 12.38.Mh, 24.85.+p 
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One of the main features of QCD, the underlying theory of strong interactions, is the spontaneous breaking of its 
approximate SU(2)xSU(2) chiral symmetry. Spontaneous breaking of this approximate symmetry explains the very 
small pion masses; in the limit of exact chiral symmetry these particles would be massless Goldstone bosons. Another 
manifestation of this is the presence of a non-vanishing quark condensate in the vacuum. It is believed that at very 
high temperature a quark-gluon plasma is formed in which chiral symmetry is restored, and much effort is being made 
O i| to explore such a phase transition by means of high-energy hadron or heavy-ion collisions. 

To describe aspects of QCD related to this symmetry it is convenient to introduce a chiral four-vector of fields 
(a, 7r), where a represents the quark condensate and the three Hi are the pion fields. In the physical vacuum, (<7,7r) 
points in the a direction. If the chiral symmetry were exact then there would be a "chiral circle" of states degenerate 
with this vacuum states. In practice the symmetry is explicitly broken by the current quark masses and so there is a 
rS \ unique vacuum. 

Because of this circle of nearly degenerate field configurations, as the chirally restored plasma cools and returns 
to the normal phase the system could form regions in which the chiral fields are misaligned, that is chirally rotated 
from their usual orientation along the a direction. There has been much recent interest in this phenomenon, which 
is known as a disoriented chiral condensate (DCC). If such a state were formed, it would lead to anomalously large 
event- by-event fluctuations in the ratio of charged to neutral pions. Since the emergence of this idea , a lot of 
theoretical work has gone into developing suitable methods for modelling the phenomenon and in exploring how it 
could be used as a signal of a phase transition in high-energy nucleus-nucleus or hadron collisions; reviews on the 
subject can be found in Refs. jl|5). 

A region of DCC can also be thought of as a coherent state of low-momentum pions. In order to produce such a 
state the hot plasma must evolve far from equilibrium and in particular it must reach an unstable configuration in 
which long- wavelength pion modes grow exponentially However it is not clear that the chiral phase transition 
is the only mechanism that could produce such a coherent pion excitation. Moreover if there were others ways of 
generating such a state, then its characteristic distribution of pions could not be regarded as a signal for formation of a 
quark-gluon plasma. In the present work we examine whether such states can form during the cooling and expansion 
of hot droplets of quark plasma coupled to chiral fields. 

Since a prerequisite for DCC formation is that the system should evolve far from equilibrium, questions of whether 
a DCC forms and how it evolves cannot be addressed in the framework of equilibrium thermodynamics but require 
either a transport theory or a hydrodynamical approach. Techniques for applying QCD directly to such situations do 
not exist at present. Hence most authors who have investigated this subject have worked within the framework of the 
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linear sigma model [0] . Some of them have made approximations that allow them to construct analytical solutions . 
Others have solved the evolution equations for the chiral fields numerically, using a variety of approaches and initial 
conditions Of these works on DCC's, the only one that, like ours, includes explicit quark degrees of freedom 

is that of Csernai and Mishustin Jl2| . Those authors studied infinite quark matter undergoing a boost-invariant 
expansion in one direction. This allowed them to look for the onset of instability with respect to fluctuations of the 
chiral fields. In contrast we study here the full evolution starting from finite-sized droplets of quark matter. 

We work within the framework of the linear sigma model ||. This is commonly used as a model for the QCD 
phase transition because it respects the SU(2)xSU(2) chiral symmetry of QCD with two light flavours of quark and 
it contains a scalar field (a) that has the same chiral properties as the quark condensate. The a field can thus be 
used to represent the order parameter for the chiral phase transition. Moreover it has been argued that the phase 
transition of the linear sigma model (without quarks) should lie in the same universality class as that of QCD 
(see also g|). 

Nonetheless universality does not provide a sufficient justification for the linear sigma model at temperatures that 
are well below the phase transition or in situations that are far from thermal equilibrium. The use of the model to 
discuss the possibility of forming a DCC thus remains an act of faith, based on the fact that the model possesses chiral 
symmetry and contains the most important degrees of freedom at low energies. If it is extended to include quarks 
(and if necessary gluons), the model can be used to describe the degrees of freedom relevant at high temperatures. At 
intermediate temperatures the quarks can also provide a rough model for the effects of massive hadrons (such as vector 
mesons and nucleons). An alternative way to include the effects of these particles would be to add terms involving 
higher derivatives of the chiral fields to the Lagrangian, as is done in chiral perturbation theory |^|. However we do 
not expect such contact interactions to provide a good description of the effects of heavy hadrons at the energies of 
interest for the question of DCC formation and we believe that inclusion of explicit quarks provides more appropriate 
way to extend the model in this context. 

A final comment that we need to make concerns the fact that the quarks in our model are unconfined and can 
escape to infinity, whereas in a more realistic model they should be converted into hadrons. We assume here that the 
hadronisation of the fast-moving quarks occurs well outside the original droplet and so does not affect the subsequent 
evolution of the chiral fields inside the droplet. 

Our basic approach is very similar to that of Ref. fl^| in that we assume a rapid quench that leaves the quarks 
out of thermal equilibrium with the chiral fields. The subsequent evolution of the system is then described by the 
classical Euler-Lagrange equations for the fields coupled to a transport equation for the quarks. The evolution of the 
quark density is described by the relativistic transport equation for fermions in the presence of chiral fields which has 
been derived by Shin and Rafelski j2^] and Zhuang and Heinz [^2] (see also |23|). Like the authors of Ref. |l2), we 
work in the classical limit where this equation has the form of a relativistic Vlasov equation for the distribution of 
the quarks in phase space. The scalar and pseudoscalar quark densities provide source terms in the Euler-Lagrange 
equations for the chiral fields. The classical system of coupled Vlasov and field equations is solved using a test-particle 
method 

In using the classical Vlasov equation we are neglecting the nonlocalities that would be present in a fully quantum 
mechanical treatment of the quarks. Our chiral fields do display strong oscillations in both space and time, but in 
general they do so in regions where the quark density is small. Hence it is unlikely that serious errors are introduced 
by the classical treatment of the quarks. We have further neglected collision terms in the transport equations that 
could give rise to dissipative effects. Although their inclusion in a transport theory for quarks poses significant 
problems p7| that go beyond the scope of our current investigation, it is important that our approach be extended to 
include collision terms since these could significantly alter the free streaming of the quarks in our present treatment. 

We study here evolution starting from a spherically symmetric droplet of quark plasma within which the quarks 
are taken to have a uniform density and a thermal momentum distribution at some temperature. The a field inside 
the droplet is taken to have the value determined self-consistently from the scalar density of these quarks. We add a 
small pionic perturbation to this configuration and then examine whether the pion field grows as the system evolves. 
Since our initial conditions do not include a realistic distribution of thermal pions we cannot estimate what ratio of 
coherent to incoherent pions is produced by the decay of such a droplet. Our purpose here is to study whether such 
a system can generate coherent pions, and if so, by what mechanism. 

Our results show that the finite size of the system plays a crucial role in determining its behaviour. In our numerical 
calculations the quarks rapidly escape from the original region in which they were placed (c/. |2^ , |2^| ). This leaves 
the chiral fields in an unstable configuration. We find that the fields always "roll" away from this configuration in 
the a direction, that is towards the true vacuum which surrounds the droplet. This demonstrates the importance of 
surface effects. In contrast to the picture suggested on the basis of bulk systems, there is no indication that such a 
droplet can form a region of more-or-less constant fields corresponding to a misaligned vacuum. Nonetheless we do 
find that, for large droplets, there can be a coherent amplification of any initial pion field. These pion fields display 
strong oscillations in space and time and so the behaviour we find is closer to that in Refs. than that proposed 
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in Ref. §. 

The mechanism that produces this amplification involves the strong oscillations of the a field that pump energy 
into oscillations of the pion field. Similar behaviour has also been seen in Refs. Such a system is thus more 

like a "pion laser" than the traditional picture of a DCC. A further difference from a DCC is that a chiral phase 
transition is not required to produce these enhanced pion fields. 

The paper is organized as follows. In Sec. |l| we present the linear sigma model that we use together with the 
classical transport and field equations. The method of solution and initial conditions are described in Sees. Ill and IV 
respectively. Our results are described in Sec. Their implications are discussed in Sec. VI and ways in which our 
approach could be improved are outlined. 



II. THEORETICAL FRAMEWORK 



A. Linear sigma model 

A simple model that embodies the main features of QCD associated with chiral symmetry is the linear sigma 
model ||. This model respects partial conservation of the axial current (PCAC) and includes a scalar field that can 
model changes in the quark condensate corresponding to restoration of chiral symmetry. The Lagrangian for the 
two-flavour version of this model is 

C = - g{<7 + in ■ t 75 )]V> + ^{d^crd^cr + d^nd^n) - U(a, n), (1) 

where a and 7r are the scalar meson and pion fields mentioned above. Although colour plays no dynamical role in our 
calculations, the quark field ip that we use describes quarks that come in three colours as well as two flavours. 
The interactions among the meson fields are described by the potential U, which we take to be of the form 

U(a,n) = ^(a*+n 2 -^) 2 -U m l*, (2) 

where /„• = 93 MeV is the pion decay constant and the parameter v 2 is given by v 2 = f 2 — m 2 /X 2 . Apart from the 
final term in U the Lagrangian is symmetric under SU(2)xSU(2) chiral symmetry. The "Mexican-hat" form of this 
potential leads to spontaneous breaking of this symmetry. In the vacuum the scalar field has a nonzero expectation 
value a = fir, which corresponds to the quark condensate of the QCD vacuum. The pions would be massless Goldstone 
bosons if the final term in U were not present. In contrast the scalar mesons have a large mass, related to the coupling 
A 2 by 

ml = 2X 2 f 2 + ml . (3) 

The mass of the a meson is often taken to be around 500 MeV, since the attractive force between nucleons can 
be described by exchange of a a meson with such a mass. However one should remember that there are important 
contributions from two-pion exchange in this channel and hence the mass of this "effective" a meson should not be 
interpreted as the mass of the underlying qq state. That state may be better identified with the /o(1370) of the 
particle data tables |50| , although Au, Morgan and Pennington |5l| have suggested that this may be part of a much 
broader structure in tttt scattering at around 1000 MeV. For most of our work we have taken m a — 1000 MeV, but 
we have also considered values in the range 600-1000 MeV. 

The quark-meson coupling constant g is more conveniently specified in terms of the dynamical mass of the quarks 
in the vacuum with spontaneously broken chiral symmetry, 

M q = gU. (4) 

The values for M q that we use lie between 300 and 500 MeV, corresponding to 3.23 < g < 5.38. 

In our studies we consider configurations in which only one component of the pion field is nonvanishing. As a 
further simplification in our numerical work we neglect the isospin dependence of the quark-pion coupling. In fact 
this coupling is unimportant for the systems we consider because the quark density is close to zero by the time that 
the pion field becomes important. In practice therefore we are working with the simplified Lagrangian 

L = - g{a + i7r 75 )]^ + \{d^a + d^d'ir) - U(a, tt), (5) 

in which the meson fields of an 0(2) linear sigma model are coupled to two flavours of quark. The extension of our 
approach to the 0(4) case is straightforward. 
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B. Equations of motion 



As in the treatment of Csernai and Mishustin ]12j , we assume that the expanding droplet of quark matter undergoes 
a rapid quench which leaves the quarks out of equilibrium with the chiral fields. We treat the subsequent evolution 
of the system classically, ignoring possible effects of quantum fluctuations. 

The classical equations of motion for the meson fields can be derived straightforwardly from the Lagrangian (Q). 
These nonlinear Klein-Gordon equations take the forms 

d^a = — A 2 (er 2 + ir 2 - v 2 )o + /,m 2 - g{^), (6) 



-A 2 (cr 2 + 7r 2 — v 2 )n — g(ipi~f5ip). 



(7) 



The quarks are described using a relativistic transport theory for fermions in the presence of scalar and pseudoscalar 
fields fin]-^. We work in the classical limit (zeroth order in K) where the distributions of quarks and antiquarks 
in phase space satisfy equations of Vlasov form. The evolution of the quark distribution in phase space f(t, r, p) is 
determined by the equation 
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V r - V r £(t,r,p) 



/(t,r,p)=0, 



where E(x, p) is the energy of a relativistic quark at the space-time point x = (t, r), 



E(x,p) = Vp 2 + M2(a;), 
and M(x) is its mass which is related to the meson fields by 



M(x) = gy/a 2 (x) +tt 2 (x). 



(8) 



(9) 



(10) 



The equation for the antiquark distribution, denoted by f(t, r, p), is identical in form to (||) since no vector fields are 
present in our model. 

The Vlasov equation (|q) can be obtained as the classical limit of an equation for the equal-time Wigner function 
of the quarks |32],[33 ' 



In fact it describes freely streaming classical quarks and antiquarks, each obeying the 
relativistic single-particle equations of motion: 
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P(t) 



E(t,r(t),p(t) 



(11) 



p(t) = -V r £?(i,r(t) J p(t)), (12) 

where the dots denote time derivatives and the dependence of the particle's energy on its position and momentum 
is given by fl|Io|). Instead of solving (g) as a partial differential equation in seven dimensions, one can replace the 
smooth distributions f(t, r, p) and /(t,r,p) by a set of classical particles obeying the equations of motion ( pT|^2| ). 
This is the basis for the test-particle method ]24|,^5| which we use to construct approximate numerical solutions to 
(||), as described in the next section. 

The couplings of a classical quark to the a and pion fields can be obtained by differentiating its energy (^) with 
respect to each of those fields. The resulting scalar and pseudoscalar quark densities are related to the quark and 
antiquark distributions by 

nil w i \ f a3 f( x 'P) + f( x >P) / 1Q % 
WW) =g<r{x) J d p — ^ , (13) 



/?• if \\ I \ f A3 /(^P) + /( a; )P) , AA s 

(i>ij 5 i>(x)) = gn(x) J d J p ' " . (14) 

These can also be obtained from the classical limit of the Wigner function, as in Refs. ||l],|2^]. Note that both scalar 
and pseudoscalar densities vanish for massless quarks (a — tt = 0). Using these expressions in the source terms for 
the fields, the field equations (|]|7]) can be rewritten in the form 
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\ 2 (a 2 (t,v)+n 2 (t,v)-v 2 )+g 2 



d 3 , 



/(t,r,p) + /(t,r,p) 



(15) 



^n(t,r) = - \ 2 (a 2 {t,v)+^ 2 {t,v)- V 2 )+g 2 J d 3 p 



(|r-v 2 MM) = - 



f{t,v, p) + /(t, r, p) 



B(t,r,p) 



7r(t,r) 



The evolution of our system is obtained by solving self-consistently the set of differential equations (g|Lf 
to an appropriate set of initial conditions (to be described in Sec. IV). 



(16) 
subject 



III. METHOD OF SOLUTION 



Solutions of the set of coupled equations (| 15 1^) can be obtained only numerically. The methods we adopt are 
very similar to those employed in Ref. |26f| to solve the similar set of equations that arise in a soliton bag model. 

The Vlasov equation (||) is conveniently solved using the test-particle method which has been widely used in the 
context of intermediate-energy heavy-ion collisions |p4| , f25f , NN scattering |^(| and evolution of a quark plasma [^8) . 
The method consists of replacing the quark and antiquark distributions by a swarm of test particles. The smooth 
distributions /(i,r,p) and /(i,r,p) are thus approximated by 



f(t, r, p) = wJ2 S a (r - r n (t)) S 3 (p - p n (t)) 



(17) 



JY 



f(t, r,p)^^^(r- f B (t)) <5 3 (p - p„(*)) 



(18) 



Each test particle follows a classical trajectory r n (t), p n (t) determined by the relativistic single-particle equations of 
motion 



r„(t) = 



Pn(t) 



E(t,r n (t),p n (t) 



(19) 



Pn(t) = -V r E(t, r n (t), p„(i)) . (20) 

Provided enough test particles are used, their distributions provide very good approximations to the smooth quark 
and antiquark distributions described by (||). In general the number of test particles is very much larger than the 
actual numbers of quarks and antiquarks we wish to describe. To account for this, the distributions must be multiplied 
by a normalization constant w. If the actual numbers of quarks and antiquarks are denoted by A and A respectively 
then the numbers of test particles are related to these by 

A= [ d 3 rd 3 p/(t,r,p) =wN, (21) 



.4 = 



J d 3 rd 3 p/(t,r,p) = wN. 



(22) 



The nonlinear equations for the meson fields (|15|,[16|) are second-order in both space and time derivatives. To solve 
these we work on a discrete mesh of points in space. In the present work we confine our attention to spherically 
symmetric initial conditions. We further require that the symmetry is maintained during the evolution, by angle- 
averaging the source terms in the field equations. In this case, the field equations can be written 



r dr 2 



(ra) - [A 2 (a 2 {t, r) + TT 2 (t, r) - ^ 2 ) + g 2 S q (t, r)] a(t, r) + / x m 2 



(23) 
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J^W- [A 2 (a 2 (t,r)+^ 2 (t,r)-^ 2 ) + 5 2 S q (i,r)] 7r(t,r), 
where the angle-averaged source term is 



s q (t, 



w 
4tt 



N 

E 



1 



-<y(r-r„(t)) 



* 1 

T — 



■*(r-f„(t)) 



(24) 



(25) 



The restriction to spherically symmetric configurations considerably lessens the computational effort required to solve 
these equations. In addition the angle average in (25|) reduces the number of test particles required to accurately map 
out the solution of the Vlasov equation. As a check on the consistency of this restriction we have run some small-scale 
three-dimensional simulations and have found no evidence for instability against non-spherical perturbations. 

The field equations are solved on a radial mesh with spacing Ar = 0.05 fm. In most of our runs we have used a 
mesh of 280 points, corresponding to a region of radius R — 14 fm. At the outer edge of the lattice we impose the 
boundary condition that the pion field vanishes while the a field is set equal to f n . This means that we do not lose 
energy from these fields when the radiated meson waves reach the edge, which is useful since we wish to calculated 
the energy radiated in particular modes of the pion field. However the fact that meson waves are reflected by the 
boundary means that our results become unphysical once those reflected waves return to the central region. For an 
initial droplet of size tq = 4 fm and a lattice of size R = 14 fm this occurs after about 18-20 fm/c. To describe 
evolution over longer times we need to use lattices with larger radii or to work with absorbing boundary conditions. 

The discrete nature of the source term due to the test particles requires that the (5-functions in ( p5| ) be smeared 
out in some way. In most of our calculations we have used a parabolic form, 



6(r- 



3 
Id 3 " 



(r-r„) 2 ), 



< d. 



(26) 



In order that all test particles contribute to the source, the thickness d should be greater than half the mesh size: 
d > Ar/2. In practice d should be somewhat larger than this to ensure a reasonably smooth source. The results 
shown in this paper correspond to d = 2Ar. We have checked that our results are not sensitive to the precise value 
of d in this region, nor do they change if Eq. (^6|) is replaced by a Gaussian smearing function. 

For the initial conditions that we study we use an initial density of test particles of about 400-500 fm -3 . For a 
region of radius 4 fm this requires a total of about 125000 test particles. We have checked that our results are stable 
with respect to a further increase in the number of test particles. 

The coupled equations are integrated using a leapfrog method with a time step of At = 0.01 fm/c. We have 
also compared our method with the more complicated staggered leapfrog method |m| described in the Appendix 
of Ref. Eq), and checked that the two agree. Since the coupled equations corresponds to a conservative system, 
conservation of the total energy provides a good test of the accuracy of our numerical algorithm. We find that the 
energy changes by less than 1% over a time of 40 fm/c (4000 time steps). 



IV. INITIAL CONDITIONS 



Our initial configuration consists of a spherical droplet of hot quark matter. We take the quark density and 
temperature to be uniform within the droplet of radius tq. The initial distribution of the quarks in phase space is 
thus 

/(0, r, p) = ^^e(r - r) r] (Eq (p) , T, M ) (27) 

where N c = 3 and Nf = 2 are the numbers of colours and flavours, respectively. 9 denotes the step function and 
T}{E,T,ii) is the Fermi distribution 

r>(E,T,»)= i + ei 1 E _ ii)/T . (28) 

For the antiquark distribution, 77 should be replaced by 

fj(E,T,n)=r)(E,T,-»). (29) 
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The quark energy appearing in these distributions has the form -Eo(p) = ^/p 2 + M§ where the quark mass Mq 
is obtained by solving self-consistently the equations for the constant a field in the presence of the (anti-)quark 
distributions ( p8| , p9| ). For a uniform quark density the equation @ for the a field takes the form 



,2N c N f 
(2tt) 3 



3 *n(E Q ,T, fj,) 



fj(E ,T, (j,) 



En 



era - Um n = 0, 



(30) 



where Eq(jj) implicitly depends on ao since Mq = gaQ. In the chiral limit (m 2 = 0) this equation always has a trivial 
solution do = 0, corresponding to a vacuum in which chiral symmetry is restored. For high enough temperatures or 
densities this is the only solution to (|30|). At low temperature and density there is also a nontrivial solution with 
lower energy which describes a vacuum in which chiral symmetry is spontaneously broken. 

The temperature at which the solution with ao =/= becomes the ground state (for fixed /i) acts like a critical 
temperature. For the parameter set m a — 1000 MeV, M q = 300 MeV this temperature is To ~ 235 MeV, while for 
m a = 600 MeV we find To ~ 150 MeV. Of course, one should remember that thermal fluctuations of the meson fields 
have been ignored in our treatment and so To is somewhat higher than the true critical temperature for the model. 
Nonetheless we can make use of the existence of this temperature to explore initial states in which chiral symmetry 
is restored as well as ones in which it is spontaneously broken. If chiral symmetry is explicitly broken then there is 
no phase transition. Instead one has a smooth crossover from states with ctq close to / w to ones with small, but still 
nonzero, values of <7q. Nonetheless, for a symmetry breaking strength that gives a realistic pion mass, this crossover 
behaviour occurs rapidly and one can still usefully define a critical temperature above which the approximate chiral 
symmetry is restored and the quarks are nearly massless. 

For the initial a field, we pick a form that interpolates smoothly between ctq inside the droplet and / T outside: 



a(t = 0,r) = o-q + (U - cr )@ a (r - r ), 



where Q is a suitably smoothed step function. In this work we have chosen to take 

6 Q (r 



ro) 



{e ar - l)/2(e ar » - 1) r<r 
(1 + tanha(r — r ) )/2 r > r , 



(31) 



(32) 



where a is determined by requiring the radial derivative of ® a to be continuous at r = ro. This is a convenient 
parametrisation whose shape is close to the self-consistent solution of the field equation for a in the presence of a 
static spherical quark distribution of the form ( ^7j ) . The constant a which describes the inverse of the surface thickness 
should be proportional to m a . 

We take the initial pion field to be of the form 



?r(* = 0,r) =7T (l-6 Q (r-r )), 



(33) 



where itq is a small, arbitrary amplitude corresponding to a fluctuation away from the self-consistent solution of (|3£ 
A nonzero initial value for cither the pion field or 7r(0, r) is needed since otherwise our system will never develop a pion 
field, as can be seen from the fact that all the terms in the field equation (|24|) contain a factor of ir(t, r). Ultimately we 
would hope to take the initial fluctuations in both a and pion fields from a thermal distribution, using an extension 
of the method described in p5[ | . However for our current investigations we use the simple ansatz ( |3l] , p3| ) in order to 
examine whether a system of this type can lead to a coherent pion field. 

The equations of motion of the meson fields (p3| , p4|) are of second order in time derivatives and so we need to specify 
the initial field velocities as well as the fields. In most of our calculations we have used static initial configurations, 
<r(0,r) = 7r(0,r) = 0. However we have also considered more general cases, and in particular ones corresponding to 
droplets that are expanding with velocity v: 



&(0,r) = -v^(0,r), 



(34) 



7r(0,r) 



(35) 
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V. RESULTS 



In this section, we present the results of our simulations corresponding to various initial conditions of the system. 
Most of our results are for our "standard" parameter set, m a = 1000 MeV, = 139 MeV and M q — 300 MeV. For 
these er and quark masses, in the chiral limit, chiral symmetry restoration occurs at the temperature To ~ 235 MeV, 
as discussed above. 

For these parameters we have examined the evolution of droplets of various sizes, with an initial temperature of 250 
MeV (above To) and zero chemical potential (zero baryon density). The initial energy density in this case is about 3.5 
GeV/fm -3 , comparable to the energy densities that are expected to be reached in ultrarelativistic heavy-ion collisions. 
This energy density is dominantly due to the quarks, the chiral fields providing only about 0.14 GeV/fm -3 . Inside the 
droplets the initial value of the sigma field, determined from the gap equation (|30|) with explicit symmetry breaking, 
is Co = 10 MeV. In order to study whether initial pionic fields can be amplified during the evolution we add a pionic 
perturbation to this configuration of magnitude ttq — 0.05/^. 

We examine first the evolution of a small initial droplet of radius ro = 1.1 fm, plotting in Fig. 1 the a and pion 
fields at several radii as a function of time. This shows the fields rolling away from the initial configuration towards 
the physical vacuum surrounding the initial droplet. Both fields then oscillate around their vacuum values, their 
amplitudes decaying as energy is carried away (mainly by long- wavelength sigma waves). Fig. 1 shows no indication 
of any enhancement of the pion field. We find similar results for small droplets with r < 2 fm. 

In contrast, for larger initial droplets we find significant amplification of the pion field, with strong oscillations in 
space and time. This is shown in Fig. 2 by the time evolution of the a and pion fields in the case of a droplet of 
radius ro = 4.0 fm, with otherwise the same initial conditions as in Fig. 1. (The fluctuations that can be seen in the 
fields at r — 0.1 fm for t < 4 fm/c are numerical noise, due to the fact that the radial density of test quarks is much 
smaller at small radii.) Another view of the same system is provided in Fig. 3, where we show the a and pion fields as 
a function of r at successive times. A larger initial volume means that the energy stored initially in the meson fields 
is larger, and this initial energy is crucial to determining whether or not the pion field is amplified. 

Also from Fig. 2 we see that the strong oscillations of the fields start only after 4-6 fm/c. This is the time needed 
by the quarks to escape the initial region. Indeed we find that for our standard choice of parameters almost all 
the quarks rapidly stream out of the initial region. This leaves the chiral fields in an unstable configuration which 
collapses towards the true vacuum. This behaviour can be seen in Fig. 4 where we plot the net density of quarks and 
antiquarks: J d 3 p (/(t, r, p) + f(t, r, p)) times r 2 as a function of r at successive times. 

In order to investigate the spectrum of the pion modes that are excited we Fourier analyse the pion field at successive 
times. From the spatial Fourier transform of the pion field, 7r(£, k), and its time derivative, ■kit, k), we calculate the 
corresponding intensity in momentum space: 

£ v (t,k) = i(|^(£,k)| 2 +c^(i,k)| 2 ), (36) 



where ujk = \Jk 2 + m 2 .. If the pion fields are sufficiently weak that they are well described by the linearised version 
of the equation of motion (|l6|), then (|3^) is simply the energy density of the pion field in momentum space. At large 
times it thus gives the energy radiated in the form of pions as a function of momentum. The total energy of the pion 
field in such a regime is 



d 3 k 



E At)= /7^3^(*> k )> ( 37 ) 



and the corresponding number of pions is given by 



^=/(0^' k) - (38) 

However, one should note that while the system is still evolving nonlinearly, ( [36|) has no such interpretation. 

Fig. 5 shows the behaviour of the intensity ([56]) at successive times as a function of momentum. (These momenta 
are the discrete values k n = mr/R, n > 1, consistent with our boundary condition at r — R.) One can see from 
this plot that during the period when the a field is undergoing violent oscillations (4-8 fm/c) the pion modes with 
momenta k ~ 0.5-2 fm -1 are significantly amplified. There is also a smaller amplification for modes with larger 
momenta, up to ~ 4 fm -1 . This enhancement of pion modes with momenta less than ~ 2 fm -1 is very similar to the 



behaviour found by Randrup 17 1, in spite of the use of a different approach. 

In Fig. 6 we plot £ n {t, k) for the most strongly amplified modes as a function of time. This shows that very little 
further energy is radiated in these modes after 8 fm/c. Similarly the total number of pions, calculated using (pq) and 
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plotted in Fig. 7, shows very little change after this time. In this case the number of pions produced is 20 times larger 
than the number corresponding to the initial pionic perturbation. 

Another related way to study the pionic modes is provided by the Fourier transform of the correlation function, 
whose time dependence was plotted by the authors of Refs. [[|||]. This quantity is equal to |7r(i,k)| 2 and so can be 
thought of as the second term of ( |36| ) divided by ui\. Once the pion field is in the regime where it satisfies a linear 
equation and oscillates harmonically 1 7r (t, k)| 2 oscillates with frequency 2cjfc. These oscillations are smoothed out 
in ( |36| ) by the presence of the time-derivative term. Although we do not present results for the correlation function 
here, we have looked at it to check that the enhanced modes shown in Fig. 6 do indeed show the expected oscillations 
for radiated pions which are in the harmonic regime. The small remaining oscillations of the pion field in the central 
region that can be seen at times after 8 fm/c in Fig. 2 do not contribute significantly to the pion radiation. This can 
also be seen from the behaviour of the total number of pions in Fig. 7. 

To explore the robustness of this mechanism for generating enhanced pion fields we have considered other initial 
conditions. These include examples with nonzero baryon density (/i ^ 0) or expanding initial droplets (v ^ in 
Eqs. (|34||35|)) and initial pionic perturbations with different magnitudes or spatial extents. In all these cases we 
obtained results that are qualitatively similar to those shown here. They are also qualitatively the same if we work in 
the chiral limit with massless pions. Our results are however sensitive to other parameters of the model. In particular 
smaller a masses correspond to smaller initial energy densities for the chiral fields. For example, with m a — 600 MeV 
a droplet of radius 4 fm leads to much less amplification than the case shown in Figs. 2-7. The effect of taking a 
larger quark-meson coupling is discussed below. 

All of the cases discussed so far correspond to evolution starting from the phase with restored (approximate) chiral 
symmetry. One can also ask about evolution from an initial state where chiral symmetry is already spontaneously 
broken. For this purpose we have considered initial conditions where the temperature is T = 200 MeV (< T ) but 
with other parameters as above. In this case, the initial hot quarks are "dressed" and have a dynamical mass of 
gaQ = 206 MeV where Co is the nontrivial solution of (|30|). We take the radius of the initial droplet and the initial 
pionic perturbation to be the same as in the previous example. The time evolution of the fields at several radii are 
shown in Fig. 8 for this case. The qualitative behaviour is similar to that of the example with T > To shown in Fig. 2. 
Again we see that a field oscillates strongly although for a rather shorter period: in this case 4-6 fm/c. During these 
oscillations there is a significant amplification of the pion field. This shows that there is nothing about this mechanism 
for generating a coherent pion field that requires a phase transition. What is essential is that the initial chiral fields 
have sufficient energy. 

Csernai and Mishustin have used the same model to study the simpler case of infinite quark matter undergoing 
a boost-invariant expansion in one direction. They found a similar coherent amplification of the pion field. In addition 
they pointed out that a more complete treatment of the model would be expected to lead to formation of clusters of 
quarks and antiquarks surrounded by domains of the coherent chiral fields. In the cases discussed above we did not 
find any such clusters (c/. Fig. 4). However if we increase the strength of the quark-meson coupling we do obtain such 
behaviour. This can be illustrated by the parameter set M q = 500 MeV (g — 5.38) and m a — 1 GeV, in the chiral 
limit. In this case symmetry restoration occurs at To ~ 170 MeV. We consider a droplet of radius ro = 2.75 fm, with 
initial temperature T = 250 MeV and chemical potential [i = 100 MeV. In Fig. 9 we show the radial dependence of 
the density of quarks and antiquarks times r 2 at successive times. In contrast to the case shown in Fig. 4, the quarks 
with smaller momenta are reflected from the surface region at around 3 fm and remain trapped within the droplet. 

The behaviour of the chiral fields shown in Fig. 10 is quite different from the previous examples. The er field in the 
central region continues to oscillate wildly around values that are much less than /„■. The quark cluster is obviously 
created in a highly excited configuration and as it settles down a significant fraction of its excitation energy is radiated 
in the form of pions, as can be seen from the fact that very large amplitude oscillations of the pion field continue for 
times up to 20 fm/c (and well beyond). In this particular example the quarks are left in a shell- like distribution with 
a radius ^2.5 fm, which is basically a "SLAC bag" soliton |56|. Other initial conditions that we have studied lead to 
more uniform distributions of quarks trapped inside shallow "bags" in the a field. In all cases the resulting excited 
quark clusters form efficient radiators of coherent pions. 

VI. SUMMARY AND DISCUSSION 

We have investigated the expansion of hot droplets of quark plasma coupled to the sigma and the pion fields of 
the linear sigma model, addressing in particular the question of whether such droplets can lead to the production 
of a coherent classical pion field or a disoriented chiral condensate (DCC). It has been suggested that the formation 
of a DCC could be used a signal for the chiral phase transition in high-energy nucleus- nucleus or hadron collisions 
. We have therefore examined whether the phase transition is the only mechanism that could produce a coherent 
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pion field. 

We treat the evolution of the system classically and solve numerically the Vlasov equation for the quarks coupled to 
the nonlinear Klein-Gordon equations for the chiral fields. These equations have previously been studied by Csernai 
and Mishustin [l2[ , but only for the case of infinite quark matter undergoing a boost-invariant expansion. Our results 
show that finite-size effects play an important role in the evolution of such systems. 

Our starting point is a spherical droplet containing a uniform density of quarks and antiquarks in thermal equilibrium 
at some temperature T. Inside this region the sigma field is taken to be constant. In using this type of initial 
configuration we are assuming a rapid quench that leaves the quarks out of thermal equilibrium with the meson 
fields. We add to this a small pionic perturbation, to act as a "seed" for possible formation of a coherent pion field. 
The fluctuations of the initial meson fields away from their self-consistent values ought to be taken from a thermal 
distribution. We plan to do this in future work using an extension of the method proposed by Randrup pq |. In 
the present paper we have used a uniform initial perturbation in order to study whether such droplets can lead to 
amplification of a coherent pion field. 

We find that the quarks rapidly escape from the original region in which they were placed, as in Ref . |28|,^9| . This 
leaves the chiral fields in an unstable configuration. Starting from the surface of the droplet, the chiral fields relax 
towards the physical vacuum which surrounds the droplet (a = = 0). The behaviour is rather different from 

the picture which has been suggested on the basis of a uniform system, namely, that the symmetry of the potential 
could allow different regions to relax towards vacuum configurations with different chiral orientations. In the finite 
systems that we have studied, the surface effects ensure that the fields always "roll" away from their unstable initial 
configuration towards the physical vacuum. 

The departure of the quarks leaves the a field in a configuration with significant potential energy. As a result it 
undergoes rather violent oscillations, during which this energy is radiated in the form of meson waves, before settling 
down to its vacuum value. If an initial pion field is present then some of this energy can be converted into pions. This 
leads to a coherent production of pions with a particular isospin and so would have similar observable effects to a 
DCC. However unlike the picture of a DCC consisting regions of differently oriented vacuum, the pion fields oscillate 
strongly in space and time. Hence the behaviour we find is closer to that in Refs. |]l|,^|, ^3|Jl^1 than that proposed in 
Ref. HI. For large initial droplets that start in the chirally restored phase, we find significant amplification of the 
number of pions, typically by factors of the order of 20 or more. The mechanism is robust against changes in the 
initial conditions as described in Sec. |y|. 

In some cases we find that the pion oscillations contain a strong component with half the frequency of the a 
oscillations. Similar oscillations of pion fields driven by the a field have also been seen in the homogeneous systems 
of a and pion fields (without quarks) studied by Boyanovsky et al. jl3| and Mrowczynski and Muller |l5|] . Analogous 
behaviour is also seen in models of inflationary universes where an additional scalar field is coupled to the field that 
drives the inflation |37| , |38|| . Mathematically one can describe this as a parametric resonance driven by the oscillations 
of the cr-dependent mass term in the pion field equation. More physically one can view it as a "pion laser" where 
energy is pumped into the pion field by the process a — > irir. The larger the initial energy in the a field, the more 
effectively this mechanism operates. Hence we find the greatest amplification for the largest droplets, in cases with 
the largest a mass. Smaller droplets, with radii less than 2 fm, do not lead to significant pion production. Also the 
enhancement factors are much smaller for cases with m a = 600 MeV compared to those where we took m a — 1000 
MeV. 

If the quark-meson coupling is taken to be strong enough, then some of the quarks remain trapped inside the droplet 
and eventually form a cold multiquark cluster, as suggested by Csernai and Mishustin (l2) . These clusters are produced 
in highly excited configurations that can act as extremely efficient radiators of coherent pions. However one should 
remember that the cold quark matter of the linear sigma model is not expected to be a very good approximation to 
hadronic matter. Hence one should not take this version of the mechanism for coherent pion production too seriously 
before more realistic models have been explored. 

The nature of the mechanism for the enhancement of the pion field found in our work makes no obvious reference 
to the chiral phase transition. As explained in Sec. |^ the model possesses a temperature Tq below which chiral 
symmetry is spontaneously broken. This enables us to study evolution from states in which the symmetry is partially 
restored, that is, where the temperature, although less than Tq, is still high enough that the a field is significantly 
reduced from its vacuum value. The matter in these states consists of dressed quarks with finite dynamical masses, 
and can be though of as a crude model for hot hadronic matter in which hadrons with higher masses than pions are 
present. We find that large enough droplets with initial temperatures below To, behave in a similar way to ones with 
initially restored chiral symmetry. In particular they also lead to amplification of coherent pions. This is a further 
indication that the mechanism involved is that of a "pion laser" which is driven by the oscillating a field and which 
does not rely on the phase transition for its operation. 

Like the DCC's and other coherent pion fields discussed in the literature, the pions produced by this mechanism 
will have the same characteristic 1/s/J distribution for the fraction / of neutral pions. However to determine whether 
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this signal will be experimentally observable over the background of incoherently produced pions we need to know 
the actual number of coherent pions. This will require the use of more realistic inital conditions. In particular, as 
already mentioned, we need to take the initial pionic fluctuations from a thermal distribution. For this purpose we 
will need to extend the method described in |35| to our model which has explicit quark degrees of freedom. In addition 
we are studying other geometries since, especially for ultrarelativistic collisions, a cylindrical system undergoing a 
boost-invariant expansion is expected to be more relevant than a spherical one. Finally we should also point out the 
need to extend our approach to include both quantum effects and collision terms in the transport equations. 
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Figure captions 

Fig. 1 : The a (a) and pion (b) fields as functions of time for various radii. The parameters of the model are: 
m CT = 1 GeV, = 139 MeV, and M q = 300 MeV. The initial conditions are: T = 250 MeV, fj, = 0, and r = 1.1 fin. 

Fig. 2 : As in Fig. 1 except for the initial radius: rg = 4 fin. 

Fig. 3 : The a (a,b) and pion (c,d) fields as functions of radius at successive times. The parameters of the model as 
well as the initial conditions are the same as in Fig. 2. 

Fig. 4 : The net density of quarks and antiquarks weighted with r 2 as a function of radius. The parameters of the 
model as well as the initial conditions are the same as in Fig. 2. 

Fig. 5 : The energy density of the pion field in momentum space as a function of momentum at successive times. 
The parameters of the model as well as the initial conditions are the same as in Fig. 2. 

Fig. 6 : The time evolution of S n for the relevant modes. The parameters of the model as well as the initial conditions 
are the same as in Fig. 2. 

Fig. 7 : The number of pions as defined in Eq. (|3^) as a function of time. The parameters of the model as well as 
the initial conditions are the same as in Fig. 2. 

Fig. 8 : The a (a) and pion (b) fields as functions of time for various radii. The parameters of the model are the 
same as in Fig. 1 while the initial conditions are: T = 200 MeV, /i = 0, ro = 4 fin. 

Fig. 9 : The net density of quarks and antiquarks weighted with r 2 as a function of radius, in the case where 
TTLfj —— 1 GeV, = 0, and M q = 500 MeV. The initial conditions are: T = 250 MeV, (i = 100, and r = 2.75 fm. 

Fig. 10 : The a (a) and pion (b) fields as functions of time for various radii. The parameters of the model as well as 
the the initial conditions are the same as in Fig. 9. 
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Abada/Birse, Fig. 1, part 1 of 2. 




Abada/Birse, Fig. 1, part 2 of 2. 
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Abada/Birse, Fig. 2, part 1 of 2. 




Abada/Birse, Fig. 2, part 2 of 2. 




Abada/Birse, Fig. 3, part 1 of 4. 
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Abada/Birse, Fig. 3, part 2 of 4. 
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Abada/Birse, Fig. 3, part 3 of 4. 
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Abada/Birse, Fig. 3, part 4 of 4. 
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Abada/Birse, Fig. 4, part 1 of 2. 
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Abada/Birse, Fig. 4, part 2 of 2. 
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Abada/Birse, Fig. 5, part 1 of 2. 
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Abada/Birse, Fig. 5, part 2 of 2. 
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Abada/Birse, Fig. 6. 
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Abada/Birse, Fig. 7. 




Abada/Birse, Fig. 8, part 1 of 2. 




Abada/Birse, Fig. 8, part 2 of 2. 




Abada/Birse, Fig. 9, part 1 of 2. 
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Abada/Birse, Fig. 9, part 2 of 2. 
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Abada/Birse, Fig. 10, part 1 of 2. 
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